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Abstract 

Double diffractive Higgs production at pp (or pp) colliders continues to attract atten- 
tion as a potential signal in the search for the boson. We present improved perturbative 
QCD estimates of the event rates for both the exclusive and inclusive double diffractive 
Higgs processes, paying particular attention to the survival probability of the rapidity 
gaps. We find that the major uncertainty is in the prediction for the survival proba- 
bility associated with soft rescattering. We show that an analogous process, the double 
diffractive production of a pair of jets with large values of Et, has an event rate which 
makes it accessible at the Tevatron. Observation of this process can therefore be used as 
a luminosity monitor for two- gluon exchange processes, such as the production of a Higgs 
boson with rapidity gaps on either side. 



1 Introduction 

A central problem in particle physics is to find a good signal with which to identify the Higgs 
boson at the present and forthcoming hadron colliders, the Tevatron and the LHC. This has 
become more important since it appears likely that the Higgs boson will be beyond the reach 
of LEP2. One possibility which, at first sight, looks attractive is to select events with a large 
rapidity gap on either side, where the conventional background is relatively low [|TJ- J5J. From 
an experimental point of view the exclusive signal looks particularly promising 

p + p — > p + gap + H + gap + p, (1) 
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and similarly for pp events. For an exclusive process we have the possibility of good experimental 
resolution on the Higgs boson mass, Mu, whereas in an inelastic collision the event rate is higher 
but the large multiplicity of secondary particles poses an additional problem in identifying the 
Higgs boson. The main question is whether the production rate of Higgs events with rapidity 
gaps is large enough. 

The cross section for double diffractive rapidity gap Higgs production can be estimated using 
perturbative QCD but unfortunately it is found || that it is strongly suppressed by rescattering 
and QCD radiative effects. Despite this, very optimistic predictions of the exclusive event rate 
persist, and are frequently cited in experimental proposals. The purpose of this paper is to 
improve the reliability of the perturbative QCD predictions by going beyond double log accuracy 
so as to give believable estimates of the event rate and to settle the present ambiguities. We will 
see that this improvement will lead to some enhancement of the event rate as compared to our 
previous estimates || [7||. Moreover, since the basic QCD mechanism for Higgs production is 
the same as that for the double diffractive central production of a pair of large Et jets, we also 
estimate the event rate for the latter process (which has a larger cross section) so that it may be 
used as a pomeron-pomeron luminosity monitor for rapidity gap Higgs production. Indeed such 
dijet data are already accessible at the Tevatron, see, for example, || and references therein. 
Thus we have a valuable check on the QCD predictions for process (0). 

As mentioned above, the cross section for Higgs production via the exclusive process (]l|) is 
suppressed by the small survival probability of the rapidity gaps. The survival probability w is 
given by the product of two factors^ 

f = S 2 spec T 2 . (2) 

First, the gaps may be filled by soft parton rescattering and, second, by QCD bremsstrahlung 
from the two fast coloured gluons which annihilate into the Higgs boson, see Fig. 1. The 
probability S^ pec not to have any extra soft rescattering was estimated in |^, [TTJ to be 

about ~ 0.1, up to a factor of 2. This suppression agrees with the simple phenomenological 
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S* pec = <e-^> * (l - — ) , (3) 

where Q is the opacity (or optical density) of the proton and (...) indicates the appropriate 
average over the impact parameter pT\ cr D is the sum of the elastic and diffractive dissociation 
cross sections and cr to t is the pp (or pp) total cross section. We will return to discuss the 
determination of S^ pec in detail in Section 5. 

The factor T 2 in (0) is the probability not to radiate gluons in the hard subprocess gg — ► H, 
and is incorporated in the perturbative QCD calculation of the exclusive amplitude. To suppress 
the QCD radiation we have to screen the colour of the annihilating gluons by an additional 
t-channel gluon, as in Fig. 1. The most optimistic scenario is to assume that this gluon, which 



1 We do not discuss the multiple (or "pile-up") interactions of high luminosities. 
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screens the colour, does not couple to the Higgs and that it has small virtuality Q\ to enhance 
the probability of screening via a large value of as- However there is a probability of relatively 
hard gluon emission coming from distance scales A > l/M^, but shorter than the characteristic 
transverse size (~ 1/Qt) at which the colour flow is screened. Such gluons fill up the rapidity 
gaps. It is found || that the typical values of Qt of the screening 'soft' gluon are indeed much 
smaller than M#, but yet are sufficiently large for perturbative QCD to be applicable. 

Since there has been much recent interest in the double diffractive Higgs signal, both for 
experiments at the Tevatron and the LHC (see, for example, [JT2[| ) , it is timely to reassess the 
estimate of the event rate. 

2 Double diffractive exclusive Higgs production 

In Ref. [|[ the production amplitude for the exclusive process ([[]), derived from perturbative 
QCD, was given by 

M( P p^p + H + p) = An* J ^e~ s (^ M ^ f( Xl ,Q T ) f(x 2 ,Q T ), (4) 

where f{x,Q T ) is the unintegrated gluon density of the proton and A is a factor associated 
with the gg —> H vertexF] 

A = (V2G F ) 1 * a s (M 2 H )/3n. (5) 

The exponential is the conventional double log Sudakov form factor which is the probability 
not to emit bremsstrahlung gluons (one of which is shown in Fig. la by px) in the interval 
Qt ^ Pt ^ Mh/2. The upper bound of pr is clear, and the lower bound occurs because there 
is destructive interference of the amplitude in which the bremsstrahlung gluon is emitted from 
a "hard" gluon k{ with that in which it is emitted from the screening gluon. That is there 
is no emission when A ~ 1/pt is larger than the separation d ~ 1/Qt of the two t-channel 
gluons in the transverse plane, since then they act effectively as a colour- singlet system. So 
the Sudakov form factor (that is the probability not to have bremsstrahlung) is given by the 
Poisson distribution exp(— S), where the mean multiplicity of bremsstrahlung is 

Sm,Ml)= £WP|)^M. (6) 

JQ\ 7T E p T 

Here E and pt are the energy and transverse momentum of an emitted gluon in the Higgs rest 
frame. 

Note that the amplitude ([|) for the exclusive process is written for forward outgoing protons, 
that is for a Higgs boson produced with small transverse momentum q?. Indeed, the presence of 
the proton form factors suppresses large qx production, and the exclusive cross section da/dyH\ 

2 We assume that Mh is below m t and hence is far from the tt threshold. 
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is calculated assuming form factors exp(6t,/2) at the proton vertices, with b = 5.5 GeV 2 . Here 
the notation means that the cross section is to be evaluated for Higgs rapidity yu = 0. 

Eq. ([|) is the perturbative QCD estimate of the exclusive amplitude to double log accuracy. 
We are now able to improve the prediction by (i) including the effect of using skewed or off- 
diagonal gluon distributions and (ii) using a more precise definition of the unintegrated gluon 
distribution. With these improvements the amplitude @ may be rewritten, to single log 
accuracy^, in the form 

M ( pp ^ p + H + p) = Avr 3 J f g ( Xl ,x[,Q 2 T ,M 2 H /A) f g (x 2 ,x' 2 ,Q 2 T ,M 2 H /A) (7) 

where f g (x, x', Q\, Mfj/4) denotes the skewed or off-diagonal unintegrated gluon density in the 
initial proton. The diagonal density is defined such that the probability to find a gluon (with 
transverse momentum Qt and momentum fraction x in the interval dQ\dx) is 
f g (dQ T /Q T )(dx/x). These unintegrated distributions are the quantities which enter when 
we apply the Qy-factorization theorem |T^| to the evaluation of the Feynman diagram of 
Fig. la. The procedure of how to calculate f g {x, x, Q T , /J, 2 ) from the conventional integrated 
gluon g(x, Q T ) is described in Ref. fH |. Here we will use the form proposed by DDT 

9 r i 
f g [x,x,Q T ,fi) = d ln g 2 [ t (Qt,V>) xg(x,Q T )\ , (8) 

where T(Qt, /i) is the survival probability that the gluon with x, x' = x and transverse momen- 
tum Qt remains untouched in the evolution up to the hard scale /z(= M#/2). T is the result 
of resumming the virtual (oc 5(1 — z)) contributions in the DGLAP evolution equation and is 
given by pi 



'HI 



dz . (9) 



The derivative dT/dlnQ 2 -, in (|8|) cancels the virtual DGLAP term in d(xg)/dlnQ' T - To be 
precise the equation for f g is a little more complicated than (§) (see eq. (3) of However in 
the relevant small x and Qt ^ Mh region, (||) is sufficiently accurate for our purposes. Note 
that after integrating (H) up to scale /i we do indeed get back the integrated gluon distribution 

J f g (x,x,Q T ^ 2 ) -^jf = T (/ i '/ i ) x 9(x,fi 2 ) = xg(x,fx 2 ). (10) 

To double log accuracy we see from (|9|) that T = exp(— S) with 5 given by (||). In fact we will 
work to single log accuracy and hence it follows from (Q) that 

T{Q T ,») = !^ML-s (11) 
ots(H) 



3 The single log accuracy of (0) may be established using (J^-factorization |13|| . The crucial point is that in 
a physical gauge (for example, the planar gauge A^n^ = with the gauge vector n M chosen as the longitudinal 
component of the Higgs 4-momentum) any additional gluon which embraces the Higgs boson (that is which 
connects the upper and lower parts of the diagram in Fig. la) gives neither a DGLAP or BFKL collincar 
logarithm. 
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We comment on this result in Section 7. 

Now we must consider the skewed effect coming from the fact that the screening gluon 
carries a much smaller momentum fraction, that is x\ <C X{. As a consequence we have 

d 



f g (x,x',Q T ,M 2 H /A) = R 9 q^ 



T(Q T ,M H /2)xg(x,Q 



(12) 



where the a/T arises because the survival probability is only relevant to the hard gluonPJ. The 
multiplicative factor R g is the ratio of the off-diagonal i'<Ci integrated distribution to the 
conventional diagonal one xg{x,Qj)- In terms of the Operator Product Expansion both the 
diagonal and off-diagonal (or skewed) distributions are given by the expectation values of the 
same conformal operators |T7|, |T8], [T9|. It was shown |2t| that for x <C 1 the expectation 
values for the diagonal and skewed distributions are the same. Hence the skewed distribution 
is completely determined in terms of the conventional diagonal gluon. Indeed for i' C i < 1 
we have HQ 

2 2A+ 3 r (A + |) 

Ra = ~W WTW (13) 

where A governs the small x behaviour of the diagonal gluon xg(x, Q^) oc x~ x . 

Note that with the VT in flT^ ) the amplitude (0) contains the same Sudakov suppression 
factor, exp(— S), as in @j. However amplitude (^) includes two improvements in comparison to 
the previous perturbative QCD form (|). We now include the dT/d\nQ T contribution, see (§), 
and also the single log part of the skewed effect, that is the factor R g in fll~2|) . Both improvements 
enhance the exclusive pp — > p + H + p event rate. This is particularly true at FNAL energies 
where the dT / d In Q\ dominates the dg/d In Q\ contribution in (§) . Moreover the enhancement 
due to R g is non- negligible. In the relevant kinematic domain we have R g ~ 1.2(1.4) leading 
to an enhancement factor Rg ~ 2(4) at LHC (Tevatron) energies. 

The values of the double diffractive exclusive Higgs cross section at Tevatron and LHC ener- 
gies that are obtained from (0) are presented in Section 4. We emphasize that this perturbative 
QCD calculation is based on the unintegrated gluon distribution f g obtained from (TT2]). It has 
been checked to be realistic in the sense that it gives reasonable cross sections for diffractive 
vector meson (p, <//?/>, Y) production |2l| and for large qt prompt photon hadroproduction at 
the Tevatron energy flUj. We will see the importance of this comment in Section 7. 



3 Double diffractive inclusive Higgs production 

The cross section for the inclusive process 

p + p -> X + gap + H + gap + F (14) 

4 It was shown explicitly in [HJ that when x' -C x only the self-energy of the hard x gluon contributes to 
the survival probability to leading log accuracy. Note that the gluon with x' ~ is almost "at rest" (that is 
Ql -C Qt) and there is no possibility of QCD radiation. 

5 Strictly speaking (|l3|) was only proved for integrated gluons However it is expected to hold equally 

well for the unintegrated distribution. 
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is much larger than for the exclusive process ([!]), see Ref. [[J. Here the initial protons may be 
broken up and so the transverse momentum of the Higgs is no longer limited by the proton 
form factor, and hence the Sudakov suppression is weaker. Unfortunately we can no longer 
achieve single log accuracy. The momenta transferred, U — (Q — fcj) 2 , are large and hence we 
cannot express the "blobs" in Fig. 1 in terms of the gluon density of the proton. At present, the 
corresponding skewed (large t«) unintegrated gluon distribution is not known. So we must use 
the (BFKL-type) non-forward amplitude. On the other hand for inclusive production we allow 
emission in a larger part of the phase space and so the QCD suppression is weaker, and the 
more approximate (essentially the double log) expression should give a satisfactory estimate of 
the event rate. 

Let us recall the main QCD formulae needed to estimate the inclusive cross section. The 
partonic quasi-elastic subprocess is ab — > a' + gap + H + gap +b r . For example, for the subprocess 
99 ~^ 9 + H + g 0, the cross section is given by 

d(J A 2 ^4 f dQ 2 dQ' 2 dtj dt 2 ( ni+n / + „ 2+n / + 5 1+5 / + 5 2+s /) /2 n ^ 

dy H A as 2% J Q2 Q ,2 h h e > ^ 

where the primes indicate quantities occurring in the complex conjugate amplitude to that 
shown in Fig. lb. Now the suppression due to QCD radiative effects comes from the double log 
resummations exp(— n.j/2) in the BFKL non-forward amplitudes, as well as from the Sudakov 
form factors exp(— S(k T , Mfj)) which arise from the requirement that there is no gluon emission 
in the interval < Pt < Mh/2. The leading logarithmic contribution again comes from the 
asymmetric configuration of the t-channel gluons, Qt fc^r. The amplitude for no gluon 
emission with Qt < Pt < i n the gap Arji is exp(— rij/2) whereQ 

--^MS)- (16) 

The inclusive cross section is obtained by convoluting the parton-parton cross sections with 
the parton densities. The results are presented in the next section. 



4 Cross sections for exclusive and inclusive Higgs pro- 
duction 

The predictions for the double diffractive Higgs production cross sections are presented in 
Table 1. The values cr cxc i given for exclusive production are obtained from the most complete 

6 There is a second (In kf T /Q^) arising from the BFKL evolution, which when resummed gives the BFKL 
non-forward amplitude <i>(A?7) exp(— rii/2). Here A77 plays the role of ln(l/a;) in the BFKL evolution. For the 
energies and rapidity gaps of interest this BFKL enhancement is small, that is <i>(Ary) « 1. Of course, a more 
precise calculation to single In Qt accuracy may give a larger amplitude due to less QCD radiative suppression. 
On the other hand the NLO ln(l/a;) corrections decrease the forward (ti = 0) BFKL amplitude. Thus we will 
still take <J>(Ar/) w 1. 
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perturbative QCD calculation that can be made at the present stage. They are based on (0) 
using the unintegrated gluon distribution of fll2|) . That is, the distributions are calculated to 
single log accuracy (as in (g) and (|])) and include the skewed effect (R g of (|TBD). We also 
include the as correction to the gg — > H vertex factor. That is A 2 given by ([5]) is multiplied 
by the regularized virtual correction p2| , or so-called i^-factor, 



\ 7T 
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1.5 A 2 . (17) 



The predictions of the inclusive cross section a- mc \ in Table 1 are obtained using the BFKL 
non-forward amplitude and are valid to double log accuracy. As mentioned above, the uninte- 
grated gluon distributions are unknown at large momentum transfer £j, and so, at present, we 
cannot improve our estimates as we have done for the exclusive case. However the values of 
<T inc i do include the factor of (|iT|). 

Finally all the cross sections in Table 1 include the survival probability S 2 pec = 0.1 arising 
from soft rescattering effects. As we shall see in the next section, the uncertainty in the value 
of S 2 pec gives by far the largest uncertainty in the predicted cross sections. 



Table 1: The cross sections a = da/dy^Q (in fb) for the central production of a Higgs boson in 
pp (or pp) collisions at y/s = 2 and 14 TeV via the exclusive or inclusive process of Fig. 1 and 
via WW fusion. The inclusive cross sections are shown for (parton level) rapidity gaps Ar] = 2, 
and also for Ar] = 3. The tabulated cross sections are obtained using S 2 pec = 0.1, but see the 
discussion in Section 5. 



M H (GeV) 


" cxcl 


^incl 


° incl (WW) 






Ar] = 2(3) 


Ar] = 2(3) 


Tevatron (y/s = 2 TeV) 








100 


0.071 


1.1 (0.09) 


0.49 (0.031) 


120 


0.030 


0.62 (0.05) 


0.41 (0.026) 


140 


0.018 


0.38 (0.03) 


0.35 (0.022) 


160 


0.008 


0.25 (0.02) 


0.30 (0.019) 


LHC (x/s = 14 TeV) 








100 


2.4 


49 (5.5) 


21.6 (5.6) 


120 


1.4 


36 (3.9) 


20.6 (5.4) 


140 


0.86 


28 (2.9) 


19.7 (5.2) 


160 


0.55 


21 (2.3) 


18.8 (5.0) 



We note that the values of cr ex ci are enhanced in comparison to our previous predictions 
which were based on the low x limit for the unintegrated gluon density. That is we took 

d{ xg{x,Q 2 T )) 
<91nQ 



assuming that \n(l/x) ^> \h{Mh/2Qt)- The predictions in Table 1, however, are based on 
the improved expression (|T2D for the skewed unintegrated gluon density. For example at LHC 
energies both logarithms (r x = ln(l/x) ~ 4.5 and ry = Ib.(Mh /2Q T ) ~ 3.5) are of comparable 
size and the derivative of a/T in (|T2"| ) implies a correction of about 1 + rx/r x ~ 1.8 to (|TB|). 
and an enhancement of the cross section der exc i oc / 4 by a factor (1.8) 4 ~ 11. Next we have 
an enhancement due to the skewed effect, R g in (|i~2"P . At first sight we might expect that the 
off-diagonal gluon distribution, 



f g (x,x') ~ yjxg(x)x'g(x'), (19) 
would be much larger than the diagonal density, f g (x, x), since x'-Ca; and x'g(x') grows rapidly 



as x' — > 0. However it was shown [^3[] that in the leading ln(l/x) limit 

f g (x,x') = f g (x,x). (20) 

Nevertheless beyond leading ln(l/x) the ratio is found to be R g ~ 1.2(1.4) at LHC (Tevatron) 
energies, leading to an enhancement R g ~ 2(4), which is included in cr exc i in Table 1. The 
third improvement is the inclusion of the single logarithmic contribution in @ for the survival 
probability T. These contributions allow for the kinematic constraints on gluon emission and 
enlarge the value of T. (It is well known that the double log expression of the type of that was 
used in 0, overestimates the suppression.) 

Table 2 shows how the QCD prediction for the exclusive cross section o" exc i would be reduced 
if we omit the various improvements one-by-one. o~ cxc i becomes 0\ if we switch off the single 
log contribution to T in @ and return to the double log formula of 0. If we then omit the 
dT I d\wQ\ term in (|| |12|) for f g the cross section reduces to 02, and finally if we omit the 
skewed effect factor R g we obtain cr 3 . 



Table 2: The reduction in the cross section cr cxc i — ► o\ — ► 02 — ► 03 due to the omission of the 
various QCD improvements one-by-one, as detailed in the text, o = da/dyH\ in fb. 





Tevatron (s/s = 2 TeV) 


LHC (v^ 


= 14 TeV) 


M H (GeV) 


120 160 


120 


160 


^cxcl 


0.030 0.008 


1.4 


0.55 


0"1 


0.012 0.003 


0.56 


0.21 


0"2 


0.8 x 10~ 4 0.8 x 10~ 5 


0.03 


0.007 


0"3 


0.2 x 10~ 4 0.2 x 10" 5 


0.012 


0.003 



We should discuss the Q\ integration which is necessary to compute the exclusive amplitude 
(0). We take the lower limit to be Q\ = 1.25 GeV 2 , since this is the starting value of the 
MRS [24] partons that we use. The saddle points of the dlnQ^ integration are at about 
<5sp — 3.2 GeV 2 and 1.5 GeV 2 for the LHC and Tevatron energies respectively. 
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The prediction for the rapidity distribution for double diffractive exclusive Higgs production 
is shown in Fig. 2. After integration over rapidity we obtainF] cr excl = 5.7 fb for Mu = 120 GeV 
at ^ = 14 TeV. 

We see from Table 1 that the values of the cross section for the inclusive process (|i~4l) are 
much larger than the exclusive cross section. Recall that the reasons are that (i) the QCD 
suppression is not so strong and (ii) the Higgs boson may be produced with a larger transverse 
momentum. However we see that for a large rapidity gap, Arj = 3, the inclusive rate falls to a 
value comparable to that for the exclusive process. 

Up to now we have dealt with a purely perturbative QCD expression. All the next-to- 
leading kinematically enhanced effects are under control and so we may hope that the higher 
order as effects will only give a 20-40% correction. Much more uncertainty comes from the soft 
non-perturbative effects. First there is the contribution coming from Q\ < Ql = 1.25 GeV 2 . 
We can estimate it using GRV partons (where we may take Qq ~ 0.6 GeV) or by freezing 
the anomalous dimension of the MRS gluon for Q T < Q , but still allowing as to run in (g). 
In both cases the contributions are comparable and not too large at the LHC energy; <7 exc i 
increases by < 20 % for M H = 120 — 160 GeV. Note that it is the inclusion of the Sudakov 
suppression effects which makes the perturbative QCD estimates infrared stable and hence the 
predictions reliable. On the other hand at the Tevatron energy the position of the saddle point 
is rather low (Q|p — 1-5 GeV 2 ) and the contribution from Qt < Qo enlarges the cross section 
by about a factor of 2. 

The main uncertainty, however, does not come from any of the above effects, but arises from 
the survival probability of the rapidity gaps with respect to the soft rescattering effects. As 
we shall see below it appears that we have been too optimistic to use in Table 1 the canonical 
value Sg. ' =0.1 at LHC energies. 



5 Suppression due to soft rescattering 

Here we will show that the main uncertainty in the calculation of the double diffractive cross 
sections arises from the estimate of how much soft rescattering fills up the rapidity gaps; that 
is, in the estimate of the probability, , for soft rescattering not to occur. 

To begin, we consider the effect of a single elastic rescattering, which is shown in Fig. 3a 
where the blob represents the whole pp (or pp) elastic amplitude including the absorptive 
corrections. It is straightforward to show that this simple rescattering amplitude gives a survival 
probability 

Sspec = i 1 ~ MBT+2b)) (21) 

7 This value may be compared to the original estimate implied by Bialas and Landshoff || of er oxc i ~ 100 fb, 
which should be multiplied by 8S^ pcc ~ 1. The factor of 8 is necessary to allow for the identity of the gluons. 
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where B e \ is the slope of the elastic differential cross section (da/dt ~ exp(B e \t)), b determines 
the t dependence of exclusive Higgs production (via the proton form factors ~ exp(6tj/2) in the 
amplitude) and cr tot is the pp (or pp) total cross section. For example at LHC energies, where 
we expect cr to t ~ 100 mb and B e \ « 20 GeV -2 , (|2"T|) gives 



S spcc « 0.65, (22) 

if we take b = 5.5 GeV~ 2 as before. 

To allow for dissociation in the rescattering process (shown by the heavier intermediate 
states in Fig. 3b) we multiply cr tot by a factor C 2 > 1, where 

cr(target dissociation) 
O = 1 ~\~ . (23) 

Here a c \ is the elastic pp cross section. Note that, for the above example, if C becomes greater 
than 1.25, then we would have already obtained a negative value for S spcc . This is a warning 
that we need to be careful in the precise values that we assume for cr to t,-Bci and C at LHC 
energies. Alternatively, we can sum up the effect of multiple rescattering using a model which 
embodies unitarity and therefore has S spec > built in. 

Let us consider the eikonal model sketched in Fig. 3c. In this model we have 

= ^ I d 2 p T (l-e-^)/ 2 ), (24) 



C 



ij / d 2 p T (l-e--(-)/ 2 ) 2 (25) 



where the impact parameter px is the transverse coordinate of the incoming proton with respect 
to the target proton, and fi(pr) is the optical density (or opacity) of the interaction. Here 
exp(— fl) reflects the absorption of the incoming beam, and exp(— CI/ 2) describes the reduction 
of the amplitude at impact parameter p?. Thus 

S spec = (e-^/ 2 ) (26) 

where the average is taken over the px dependence, exp(— p T /4b), of the amplitude for exclusive 
Higgs production 

M oc e fc( ' 1+t2)/2 . (27) 

That is we have 

J d 2 PT e-^V 2 e~^ b 

J d 2 p T e~^ b • 1 j 

For the opacity we take the Gaussian form 

%) = ^e^, (29) 
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where a(s) = <j (s/s ) a corresponds to the Pomeron exchange amplitude shown by the double 
lines in Fig. 3c. We take the slope of the Pomeron amplitude to be 

- = B + 2a' P ln(s/so)- (30) 

We tune the parameters (a ,B ,A, a' P ) of the eikonal model to describe the behaviour of a tot 
and the pp elastic scattering data throughout the ISR to Tevatron energy range (30 < yfs < 
1800 GeV). Finally we study the predictions for S spec for two relevant values of the enhancement 
parameter C of (|23f) , namely C = 1.15 and C = 1.3. The smaller value of C is obtained if 



we include only the nucleon resonance excitations [p7 |; in terms of partons it means that we 
account mainly for valence quark rescattering. On the other hand at the larger (LHC) energies 
we have to include rescattering of partons with small x ('wee' partons), and in this case C ~ 1.3 
is more appropriate. Both choices allow a satisfactory fit of cr tot and the elastic dataQ, and the 
values of S^ pcc obtained from (PS|) are shown in Table 3. 



Table 3: The probability S^ pec , that the rapidity gaps survive rescattering, calculated using the 
eikonal model fl28D for two values of the enhancement factor C of (^3|), namely C = 1.15 and 
C = 1.3 (expected to be appropriate for Tevatron and LHC energies respectively). The slope 
b for exclusive Higgs production, ([27]), is expected to be 5.5 GeV -2 , but smaller values are not 
excluded (see text). 



b GeV" 2 


Tevatron (y/s 


= 2 TeV) 


LHC (Vs 


■= 14 TeV) 


C=1.15 C = 1.3 


C= 1.15 


C= 1.3 


5.5 


0.11 


0.034 


0.054 


0.011 


4 


0.07 


0.013 


0.029 


0.003 


2.5 


0.04 


0.003 


0.012 


0.0003 



From the results shown in Table 3 we see that the survival probability S^ pcc depends sensi- 
tively on the value of the slope b (of (p7[)), that is on the spatial (p?) distributions of partons 
inside the proton, see (|28|). Unfortunately we cannot be completely sure that the value we have 
adopted, b = 5.5 GeV -2 , is correct. This value is obtained by assuming that the distribution 
of colour dipoles (gluons) mimicks the electric charge distribution of the proton. However in 
diffractive J/ ip electroproduction the slope is observed to be ~ 4 GeV -2 |28| . Since this pro- 



cess is also mediated by two-gluon exchange, another choice for b could be 4 GeV 2 . Moreover 
the J /if) slope is given by 

b^ = b + 2a' P ln{l/x il ), (31) 

where l/x^ w W 2 /M^. Thus, since the J /if) HERA data |28f sample x ~ 10 3 , whereas Higgs 
production at the LHC samples x ~ 10" 2 , it suggests that b is about (2a P lnlO) less than b^. 

8 The parameter A, which specifies the Pomeron intercept, is found to be A = 0.10 and A = 0.13 in order to 
fit the "soft" data, taking C = 1.15 and C = 1.3 respectively. 
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For this reason we also show results in Table 3 for b = 2.5 GeV~ 2 , which if it were true would 
give much lower values of S 2 ■ However the evidence from HERA is conflicting. Large Q 2 
open qq or p meson electroproduction |2S] are observed to have slope b p ~ 7 GeV -2 , in the 
same region of x as J/ip. From this point of view the lower J/^-motivated values of b look 
anomalous. So, on balance, the value b = 5.5 GeV~ 2 looks to be the most realistic. 

Taking b = 5.5 GeV -2 , we see from Table 3, that at the Tevatron we have S 2 C ~ 0.1, 
which is the canonical value that we have used. However at the LHC^] the corresponding value 
is S 2 pec ~ 0.05. Even worse, for the more relevant choice of enhancement factor, C = 1.3 at 
LHC energies, we predict S 2 pec ~ 0.01. 

We have used other models to estimate the survival probability S 2 pec and, given the values 
of C and b, we have found essentially the same results as in Table 3. The reason is interesting. 
At high energies the centre of the proton becomes black, that is Q ^> 1 and exp(— Q/2) ~ 
in Hence the main contribution to S spec comes from the peripheral region ^ Po, where 
Po is where the proton starts to become transparent, i.e. fl(pr ^ Po) ~ 1- Thus, as long as the 
model describes cr to t and elastic pp data, the prediction for S spec does not depend crucially on 
the details, but is controlled essentially by the proton radius (or B c \ « er 2 ot /167rtT e i) and the 
slope b. 

Let us finally comment on S 2 pec for a Higgs boson produced by WW or 77 fusion. It is 
not excluded that the radius of the quark distributions in the proton is larger than that of the 
gluons. If this were the case, then S 2 pec for WW fusion would be larger than those shown in 
Table 3. The most exciting example is 77 — > H production. This process takes place at very 
large impact parameter p T , and here S spec ~ 1. In Ref. the 77 — > H cross section at the 
LHC was predicted to be 0.3 fb for Mh ~ 120 GeV. This would be close to our prediction of 
0.6 fb for production by two-gluon exchange if we were to take a survival factor of S 2 C = 0.01, 
instead of S 2 pcc = 0.1. If indeed this is the case and, moreover, if we were to assume that 
b < 5.5 GeV -2 is correct, then we would have the following hierarchy 

1 « S 2 pec (77-#) » S 2 spcc (WW - H) > S 2 spcc (PP - H), (32) 

where PP denotes the two-gluon exchange mechanism of Fig. 1. Of course for the default 
choice b = 5.5 GeV~ 2 (i.e. assuming the spatial distributions of quarks and gluons to be the 
same) we have 

Sl ec (WW - H) = S 2 spcc (PP - H). (33) 



6 The dijet monitor 

The previous section demonstrates that the 'Achilles heel" of the calculation of the Higgs 
production cross sections is the uncertainty in the soft survival factor S 2 pcc . Fortunately there 

9 Note that the simple formula (||) would give S'spoc independent of collider energy if a D /crtot were constant, 
which is not incompatible with the present data for a D /atot, although the errors are large. However this 
estimate is too naive since it assumes the same spatial distribution of partons in the proton for both the soft 
and hard processes, that is b = B c \/2. 
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is a way to experimentally measure S^ pec by observing the double diffractive production of a 
pair of high E T (~ M h /2) jets with rapidity gaps on either side of the pair. The process is 
described by the same Feynman diagrams, both in the case of the exclusive process (|I]) and 
also for inclusive production (|T3). Essentially we need simply to replace the gg — ► H subprocess 



by that for gg — > dijet. The dijet event rate is much larger than that for the Higgs signal and so 
the collider experiments should be able to directly test the QCD estimates and measure S^ pec . 

QCD estimates of the rapidity gap dijet rate were given in 0. In Table 4 we present 
improved numerical results in a kinematic range comparable to that for Higgs production. We 
use the same prescription that was employed to calculate the Higgs production cross sections 
presented in Table 1. That is for exclusive dijet production we integrate from Q T = Q to 
Qt = Et over skewed unintegrated gluons (0), calculating the QCD radiative survival factor 
T to single log accuracy. However the NLO /Y-factor for the subprocess gg — > dijet is omitted. 
This correction depends on the "jet finding" algorithm. Usually the size of the jet cone is chosen 
in such a way that the effective NLO /Y-factor is close to 1. For comparison with Table 1 we 
continue to use the canonical soft rescattering factor S% =0.1, although we note from Section 
5 the true factor may be smaller. 

In practice it is impossible to study a purely exclusive dijet production process, analogous 
to (H|). We cannot distinguish a bremsstrahlung gluon emitted in the dijet rapidity interval 
from a gluon which belongs to one of the high-^ jets. We have therefore chosen rapidity 
gaps such that bremsstrahlung is only forbidden for \rj g \ > 2 in the dijet centre- of-mass frame. 
Of course the QCD radiative suppression will be much smaller in this case. For example, for 
semi-exclusive production of E T = 50 GeV jets at LHC energies, we have a typical survival 
factor T(Q SP ,M H /2) ~ 0.5 at the saddle point Qg P ~ 2 GeV 2 of the dhiQ T integration. 

The double diffractive dijet cross sections are much larger than those for Higgs production. 
For example if we take a dijet bin of size SEt = 10 GeV for each jet and 8(j]i — 772) = 1 we 
estimate, for Et = 50 GeV jets at LHC energies, 

dcr exc i/ drj\ ~ 38 pb, da mc i/dr]\ ~ 240 pb (34) 

where 77 = (771 H- 772) /2, and the rapidity gaps are taken to be A^(veto) = (r7 min , 77 max ) = (2, 4.1) 
for the inclusive case (see @ for the definition of the dijet kinematics). For 30 GeV jets at the 
Tevatron the corresponding cross sections are 

da exd /dr]\ ~ 17 pb, da ind /dr]\ ~ 150 pb. (35) 



The numbers in square brackets in Table 4 correspond to double diffractive bb dijet pro- 
duction. For Mh in the range that we consider, this process is the main background to the 
double diffractive Higgs signal. However the event rate of bb jets is more than two orders of 
magnitude lower than the gluon dijet rate. Even after integration over a 5Et = 10 GeV interval 
the rate is comparable to the Higgs cross section. We conclude the bb background should not 
be a problem. 
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We emphasize that the rapidity intervals chosen in our calculations refer to partonic rapidi- 
ties. Table 4 shows that the cross sections depend sensitively on the size of the rapidity gaps, 
Ai] vcto . Practical estimates require Monte Carlo simulations appropriate to the specific experi- 
mental cuts and which include treatment of possible initial state radiation and the hadronization 
of the large Et jets. 

Table 4: The cross sections a = da / ' dp T dr]dAr]\ (in fb/GeV 2 ) for the exclusive and inclusive 
double diffractive dijet production, for three values of the rapidity difference of the two jets 
Arj = — 772. The results are shown for different collider energies (y/s) and different transverse 
energy (Et) of the jets. The two jets are taken to have the same Et- The rapidity gaps are 
defined by the intervals ±A^ vcto = ±(r/ min , ^ max ) JTJ. The numbers in the square brackets are 
the bb component of the dijet signal. 





A77 


0"cxcl Hi) 


O-mclO'j) 


[flWi(66)] 




[cr exc i(66)] 


Ar/ veto = (2,4.1) 


Ar] vcto = (1.5,4.6) 


y/s = 2 TeV 





0.97 [0.008] 


5.2 [0.04] 


0.28 [0.002] 


E T = 50 GeV 


1 


0.76 [0.005] 


4.4 [0.03] 


0.23 [0.0015] 




2 


0.31 [0.001] 


2.4 [0.01] 


0.11 [0.0004] 


y/s = 2 TeV 





29 [0.23] 


240 [1.9] 


15 [0.12] 


E T = 30 GeV 


1 


22 [0.13] 


220 [1.5] 


13 [0.09] 




2 


11 [0.03] 


140 [0.5] 


8 [0.3] 


y/s = 14 TeV 





38 [0.30] 


240 [1.9] 


24 [0.19] 


E T = 50 GeV 


1 


31 [0.22] 


240 [1.6] 


23 [0.16] 




2 


19 [0.08] 


200 [0.7] 


18 [0.07] 



7 Comparison with other QCD-based predictions 

We have argued that perturbative QCD gives reliable estimates for double diffractive Higgs and 
dijet production, up to the uncertainty in the soft rescattering effects. It is therefore important 
to understand the origin of the difference with recent more optimistic estimates of the event 
rates. 



Double diffractive high Et dijet production has been recently estimated by Berera [30|. His 



non-factorized N(L)DPE model is similar to our perturbative QCD approach. However there 
are some differences in application. First, in |30| a fixed coupling as{E T ) was used in the 



double log form of the Sudakov suppression factor, whereas here (and in Ref. J7|) we use the 
more appropriate running coupling ols(Pt) inside the integration Q T < p T < E T . Second, dijet 
production cannot be a pure exclusive process since it is impossible to forbid extra emission 
in the central rapidity interval occupied by the dijets. Thus we have, at best, a semi-exclusive 
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reaction in which the suppression^ is only relevant in some rapidity interval Srj veto , see Ref. 0. 
When this is taken into account the predictions |30] are reduced by an amount that roughly 
compensates for the enhancement due to the use of as{Ej). The most important numerical 
difference between [30| and our predictions is due to treatment of the gluon exchanges. In 



Ref. |3(| (and also in [p5fl , as we discuss below) a non-perturbative two-gluon model is used 
in which the gluon propagator is modified so as to reproduce the total cross section. On the 
contrary, we have used a realistic unintegrated gluon density, determined from conventional 
gluons of global parton analyses, which has been found to give a consistent description of other 
processes described by perturbative QCD. 

It was emphasized in |3(| (see also [0) that the above non-perturbative normalisation based 
on the value of the elastic (or total) cross section fixes the diagonal gluon density at x ~ It I 1 \fs 
where the transverse momentum is small, namely £t < 1 GeV < Qq. Thus the value of x is 
even smaller than 

x « Qt/Vs < x « Mn/yG. (36) 

However, the gluon density grows as x — > and so it is clear that such a non- perturbative 
gluon normalisation will overestimate the double diffractive cross section. 

We now turn to the very optimistic ( "brave" ) estimate for double diffractive exclusive Higgs 
production, process (|IJ), that has recently been presented in |2"5]| . For M# = 100 GeV the 



prediction is da/dy ~ 20 fb, even for Tevatron energies, see Fig. 2. We are unable to justify 

par 



this estimate. First the radiative suppression T 2 (= S* 2 in fl25l ) ~ 0.1 which is much larger than 



our determination. The phenomenological estimate of T 2 in p5j is based on the known hadron 
multiplicity A^ a d measured in the rapidity interval Ay = ln(M#/so) in a soft hadron-hadron 
collision. This multiplicity increases as ln(M^/so) and has nothing to do with the double 
logarithmic bremsstrahlung, where the mean number of emitted gluons S oc In 2 M\. Next, 



motivated by the BLM prescription [32|, the coupling as in the gg — > H matrix element (that 



is in the analogue of ([5])) is evaluated in [25] at a low scale Qo rather than Mh- However to apply 
the BLM procedure consistently we must ascertain which part of the gluon self-energy insertions 
are already included in the survival factor T (or S 2 ), and which part should be attributed to 
as- Indeed in calculating T to single log accuracy we obtained the pre-exponential factor 



ois{Q 



as(M 2 /4)' 



(37) 



see (O). This factor reflects the fact that, as usual, the double log approximation overes- 
timates the kinematically available phase space for emission. That is the probability not to 

10 So, the Sudakov suppression appears to be much weaker for dijet production than for the analogous Higgs 
process. We can easily gain insight into the origin of this difference by recalling a similar phenomenon in 
the radiative effects accompanying the production of narrow and wide heavy resonances (see, for example, the 
discussion of well known QED effects in Ref. J3l[]). The energetic bremsstrahlung pushes the initial state off 
the resonance energy for the non-radiative process. Thus, the narrow resonance could be produced only if we 
damp radiation with energy exceeding the resonance width. The wider the resonance, the larger is the phase 
space available for emission and, therefore, the less pronounced is the Sudakov suppression. 



15 



bremsstrahlung a gluon, T, is larger than exp(— S). Thus, in conclusion, as far as the single 
log corrections are already included in the T factor, we must use as{Mfj) in (|), together with 
the K-f actor of ( |TTD evaluated at scale /i = M H . 

From Table 1 and Fig. 2 we see that the perturbative QCD predictions for da/dyn show 
a strong increase with increasing energy, which arises because of the growth of the gluon 
densities xg(x,Q^) with increasing 1/x ~ s/Mfj. On the contrary the non-perturbative two- 
gluon-exchange-type phenomenological models have no x dependence. The predictions of these 
models depend only weakly on energy through the energy dependence of the "soft" cross section 
which is used to normalise the two-gluon exchange contribution. The same arguments apply 
to the production of a pair of high jets. Therefore an experimental study of the dijet 
production rate as a function of the collider energy will clearly be able to discriminate between 
the perturbative QCD determinations and the non-perturbative model approaches. 

8 Summary 

We have calculated the cross sections for exclusive and inclusive double diffractive Higgs boson, 
and also dijet, production in the central region, at both LHC and Tevatron energies. That is 

pp — > p + (H or jj) + p 

(38) 

pp — > X + (H or jj) + Y 

where + denotes a rapidity gap. These processes are driven by 'asymmetric' two gluon exchange, 
with the colour screening gluon being comparatively soft, but still in the perturbative QCD 
domain. 

All the important perturbative QCD corrections were included in the calculation. A pre- 
scription for the unintegrated gluon distribution, up to single log accuracy, was used. The major 
uncertainty comes from the non-perturbative sector, namely from the value of the survival fac- 
tor Sg pec — the small probability to have no secondaries from soft rescattering populating the 
rapidity gaps. We found that this probability S^ pec depends sensitively on the spatial distribu- 
tion of gluons inside the proton. In the tables in which we presented the Higgs and dijet cross 
section predictions, we took S^ pec = 0.1, but at LHC energies our estimates of S^ pec indicate 
that the value is most likely to be an order of magnitude smaller. 

To overcome the normalisation uncertainty due to the lack of knowledge of we proposed 
that measurements of the double diffractive production of dijets would act as a luminosity 
monitor for the two-gluon exchange processes of The estimates of the dijet cross section 
are such that the process should be readily observable at the Tevatron and at the LHC. In 
particular measurements of jets with E T ~ M H /2 would enable the cross section for the double 
diffractive production of the Higgs boson to be reliably predicted, since the two processes are 
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driven by the same S spec factor in the same kinematic region. Unfortunately even our most 
optimistic predictions for the Higgs process are considerably smaller than previous estimates, 
and would make the process hard to observe at the Tevatron and the LHC. 
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Figure Captions 



Fig. 1 Diagrams for (a) exclusive, and (b) inclusive, double diffractive Higgs production of trans- 
verse momentum q<? in high energy pp (or pp) collisions. The QCD radiative corrections 
(such as the emission of the gluon of transverse momentum px) suppress the number of 
rapidity gap events via Sudakov form factors, exp(— S). For inclusive production qx can 
be much larger and the Sudakov suppression is weaker; however there are additional QCD 
radiative effects from the double log resummations exp(— rii/2) in the BFKL non-forward 
amplitudes. 

Fig. 2 The perturbative QCD predictions of the rapidity distribution for double diffractive exclu- 
sive Higgs production at LHC and Tevatron energies. We also show the recent prediction 
by Levin [p5] , which is discussed in Section 7. 



Fig. 3 Diagram (a) illustrates the absorptive correction to exclusive Higgs production, assum- 
ing that only elastic pp rescattering occurs, with an amplitude Iul4 = sutot- Diagram 
(b) includes both elastic and inelastic intermediate states. Diagram (c) is an eikonal 
representation of (b), where the double line denotes Pomeron exchange. 
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(a) Exclusive production (b) Inclusive production 




Fig.l 



